Emergent electromagnetism in solids 
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Electromagnetic field (EMF) is the most fundamental field in condensed-matter physics. Interac- 
tion between electrons, electron-ion interaction, and ion-ion interaction are all of the electromagnetic 
origin, while the other 3 fundamental forces, i.e., gravitational force, weak and strong interactions 
are irrelevant in the energy /length scales of condensed-matter physics. Also the physical proper- 
ties of condensed-matters such as transport, optical, magnetic and dielectric properties, are almost 
described as their electromagnetic responses. In addition to this EMF in the low energy sector, it 
often happens that the gauge fields appear as the emergent phenomenon due to the projection of 
the electronic wavefunctions onto the curved manifold of the Hilbert sub-space. These emergent 
electromagnetic fields (EEMF's) play important roles in many places in condensed-matter physics 
including the quantum Hall effect, strongly correlated electrons, and also in non-interacting electron 
systems. In this article, we describe its fundamental idea and some of the applications recently 
studied. 

PACS numbers: 73.43.Cd,72.25.-b,72.80.-r 
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I. INTRODUCTION AND BASIC CONCEPTS 



Electromagnetism is the first gauge theory which was 
recognized in physics, and essentially contains the prin- 
ciple of the relativity. Most of the physical properties 
of solids are regarded as the responses of the systems to 
the electromagnetic field (EMF) . For example, the trans- 
port phenomena are described by the currents induced 
by the electric field, magnetism is the behavior of the 
magnetization in response to the magnetic field, and op- 
tical properties are described by the electric polarization 
driven by the EMF of light. On the side of the electronic 
system, the electron-electron and/or electron-ion interac- 
tions mediated by the EMF, i.e., mostly the Coulomb in- 
teraction, determine the quantum state even without the 
external EMF as a probe. This interaction often leads to 
the emergent electromagnetic fields (EEMF's) as will be 
described in more details below. The other aspect of the 
electromagnetism in solids is the various cross correla- 
tion effect. Namely various different combinations of the 
inputs and responses have been explored. For example, 
an applied magnetic field can induce the electric polar- 
ization, or the magnetization emerges under the electric 
field in the magneto-electric (ME) effect. This effect was 
proposed long time ago [l| , but has been very weak. An 
explosion of the research was triggered by the discovery 
of the multiferroic materials showing both the magnetic 
order and ferroelectricity which are coupled strongly to 
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each other 0, Q. These phenomena are driven by the 
internal degrees of freedom of electrons, i.e., spin, charge 
and orbital, which are again described coherently by the 
EEMF's. 

This EEMF is analogous to the usual EMF, but is 
much richer because of the following aspects [3, [B[ ■ (i) 
It is defined on the lattice because of the crystal struc- 
ture and Bloch waves in solids, and the lattice gauge 
theories are relevant. Therefore, various topological de- 
fects play important roles, (ii) One can generalize the 
space where EEMF is defined including the real space, 
momentum (k-)space, and other parameter space, (iii) 
The gauge group is not only U(l) as in the case of EMF, 
but can be more general such as the SU(2) corresponding 
to the Kramer's doublet in the time-reversal symmetric 
system, and U(N) corresponding to the N-bands system, 
(iv) Topological terms such as Chern-Simons term, WZW 
term, and #-tcrm often appear to the effective action for 
EEMF. Therefore, solids can be an ideal arena to study 
various concepts and techniques in quantum field the- 
ory 6]. 

The basic reason why the EEMF's appear in 
condensed-matter is the restriction of the Hilbert space. 
Namely, only the small portion of the Hilbert space is 
relevant to the low energy phenomena, which constitutes 
the fiber-bundle structure. For example, consider a vec- 
tor moving on the surface of unit sphere S 2 . The par- 
allel transport of this vector on S 2 leads to the change 
in its direction as shown in Fig.l. This is because S 2 
has the (constant) curvature and the vector picks up this 
curvature through the connection for infinitesimal paral- 
lel transport. This connection corresponds to the vector 
potential of EMF while the curvature to the magnetic 
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flux, which are related via Stokes theorem to each other. 
An analogy can be drawn between this geometry prob- 
lem and the quantum mechanics Q. The vector cor- 
responds to the wavefunction and S 2 to the restricted 
Hilbcrt space. Therefore, the parallel transport on the 
subspace leads naturally to the gauge structure. We dis- 
cuss below briefly the Berry phase [8| and the spin-orbit 
interaction from this viewpoint. 



A. Berry phase of electrons in perfect crystal 

The Bragg diffraction of the electrons by the periodic 
electric field of the ions in a crystal leads to the Bloch 
states as described by 



^nk<r(r, S) 



(1) 



where w n k(r) is the periodic function with respect to the 
translation of lattice vectors, x<r( s ) * s the spin wavefunc- 
tion for a =1\4- with the spin coordinate s = ±. The 
band index n and the crystal momentum k are the good 
quantum numbers and the energy eigenvalues e„(k)'s 
with different 7i's are usually separated by the band gaps. 
This means that the electrons in a solid are confined to 
each band in the lower energy sector than the band gaps, 
which corresponds to the projection of the wavefunctions 
to the sub-space of the Hilbert space. This sub-space is 
usually " curved" and is represented by the gauge connec- 
tion. (This is an application of the Berry connection to 
Bloch electrons ||.) It is given by the expression 



a na (k) = -i(u n k\dk a \u nk ) 



(2) 



with a = x,y,z. This quantity appears in the inner- 
product of the 2 wavefunctions for the 2 neighboring k- 
points as 



(wnk|wnk+Ak) = cxp(ia„(k) • Ak). 



(3) 



This Berry connection leads to the concept of Berry cur- 
vature corresponding to the magnetic field in the momen- 
tum space as 



b na (k) = e abc d kb a nc {k) 



(4) 



with the totally antisymmetric tensor e a bc- 

What is the physical consequences of this Berry phase 
? To answer to this question, it is useful to consider 
the motion of a wave-packet made of the Bloch states, 
where the position and momentum of a particle is de- 
fined within the accuracy consistent with the uncer- 
tainty principle 9]. In the presence of the Berry con- 
nection, the position operator x originally defined as 
x u = id/dk u should be generalized to the gauge covari- 



ant form 



id/dkfj 



(k). This is dual to the 



momentum operator tt u in the presence of the magnetic 
field, i.e., w u = p u + eA u (x.) = -id/dx^ + eA ft (x), where 
A M (x) is the vector potential of the EMF. Analogously 
to the commutator [7r x ,7T 9 ] = iaB z (x) with -B z (x) being 



the z-component of the magnetic field, one can derive 
[x,y] — ib z (k). Namely, the real space coordinates do 
not commute with each other. With this commutator, 
the equation of motion for the wavepacket reads [§, Hfjj 



~dt 

dt 



(5) 



Here the Hamiltonian H is given by H — e„(k) + V^(x) 
with £n(k) being the energy dispersion of the band n of 
interest, and V(x) is the slowly varying external poten- 
tial. Putting this into eqs.([5]) and we obtain 



dx n 
~dT 



d£n(k) 

dk n 
de n (k) 



, 9V{yi) 



dx v 
dV{x) 



(b x F) M , 



<9e„(k) 



dk a 
~dT 



dk a 
dV(x) 
dx u 



(6) 



where F = — V x y(x) = eE is the force acting on the 
electron. In eq.©, the duality in the equation of mo- 
tion is evident, i.e., the real space quantities x a ,A u {-x), 
V(x) and those in the momentum space fc M , a na (k), e n (k) 
enter into the equations of motion in a symmetric way. 
Therefore one can translate the phenomena in real (mo- 
mentum) space to those in the momentum (real) space. 
Note, however, that compared with the real space mag- 
netic field which is divergence- free, i.e., V x • B(x) = 0, 
corresponding to the absence of the magnetic monopole 
in real space, Vk • b(k) can be non-zero in momentum 
space. Another important remark is that the symmetries 
give the following constraint. The time-reversal symme- 
try T gives the relation b n (k) = — b„(— k), while the 
inversion symmetry / b„(k) = b„(— k). Therefore, when 
both T and / symmetries are there, there is no Berry 
curvature b„(k). Also, in the noncentrosymmetric sys- 
tem where /-symmetry is absent, b„(k) can be nonzero 
although the contributions from k and — k cancel with 
each other. 



B. Emergent SU(2) gauge field from Dirac 
equation 

Starting from the Dirac equation, the projection of 
the wavefunctions onto the low energy subspece leads to 
the gauge structure analogous to the electromagnetism. 
This formalism is based on the real-space picture, and is 
complementary to the discussion in the last subsection 
where the momentum space geometry has been consid- 
ered. When the SU(2) spin space is preserved, it leads to 
the SU(2) non-Abelian gauge field which is coupled to the 
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spin current, corresponding to the spin-orbit interaction 
(SOI). 

We start with the quantum electrodynamics (QED), 
where the Dirac relativistic electrons and their charge 
current is minimally coupled to EMF. Therefore, it ap- 
pears that there is no chance for the spin current to take 
some role in the electromagnetic phenomena. More ex- 
plicitly, in the natural unit where % = c = 1, the La- 
grangian of QED reads @ 



L — iplij^Dp — m]ip, 



(7) 



where ip is the 4-component spinor field, ip = ip'^f , 7 M 
is the Dirac matrices, and — — ieA^ is the gauge 
covariant derivative with fi = 0, 1, 2, 3. The 4-component 
charge current density is defined as 



r 



dL 



(8) 



whose 0-component is the charge density p, while the 
spatial components are the current density j. From the 
gauge invariance, the conservation law of the charge is 
derived through Noether's theorem as 



V = | + v.j = o. 



(9) 



By taking the variation, one can derive the Maxwell equa- 
tion 



and Dirac equation 



[ij^Dn -m]ip = 0. 



(10) 



(11) 



As is well-known, the solutions to Dirac equation are clas- 
sified into 2 classes, i.e., the positive energy and negative 
energy states separated by the twice of the rest mass en- 
ergy of the electrons 2mc 2 . Since the energy mc 2 is of 
the order of MeV, for the low energy phenomena typi- 
cally of the order of ~ eV, the negative energy states are 
not relevant. Therefore, the non-relativistic Schroedinger 
equation is usually used which describes the dynamics of 
the 2-component spinor wavefunctions for positive en- 
ergy states. However, one needs to take into account 
one important aspect of "projection". Namely, the neg- 
ligence of the negative energy states means the projec- 
tion of the wavefunctions onto the positive energy states, 
i.e, the sub-Hilbert space. Usually the sub-space is not 
flat but curved, and associated geometrical structure is 
introduced. The derivation of the effective Lagrangian 
describing the low energy physics is achieved by the ex- 
pansion with respect to l/(mc 2 ), and the result reads 



D 2 



1 



iip^ Do^p+ip^ tj)-\ ip' 

2m 2m 



eqr a A- A a +—A a -A a 

4 



(12) 



where i\> is now the 2-component spinor and Dq = 
do + ieA + iqA^^Y, and Z?.j = di — ieAi — iqAf 1 ^- 
(i = 1,2,3) are the gauge covariant derivatives with 
q being the quantity proportional to the Bohr magne- 
ton [III is the vector potential for EMF, while 
the SU(2) gauge potential are defined as Aq = B a , 
Af = eiaiEg. The former is coupled to the charge cur- 
rent, and the latter to the 4-component spin current 
jg = i>a a ^p, jf = ^[^a a D^ - D^G a ^\. Note that 
the spin current is the tensor quantity with one suffix for 
the direction of the spin polarization while the other for 
the direction of the flow. Note an important difference 
between the EMF and the SU(2) gauge field. The former 
has the gauge symmetry, i.e., the freedom to chose the ar- 
bitrary gauge for the vector potential A^ , while the " vec- 
tor potential" A a ^ for the latter is given by the physical 
field strength B and E. Actually, the relation d^A^ = 
holds. Therefore, the SU(2) gauge symmetry is absent. 
This is the basic reason why the spin is not conserved 
in the presence of the relativistic SOI. Instead, the spin 
current is " covariantly" conserved and satisfies [n]-[l3| 



D Jg + D- J a = 0. 



(13) 



This means that in the co-moving frame the spin is con- 
served, while in the laboratory frame the spin source or 
sink appears when the electron forms a loop and comes 
back to the same position in space since the frame already 
changes [l3j]. (Note that the usual SU(2) gauge theory 
is a nonlinear theory and the gauge field is "charged", 
and the sum of the spin current by matter field and the 
gauge field is conserved in the non-Abelian gauge theory 
as Yang-Mills first showed [6|. However, the SU(2) gauge 
invariance is absent here, and also the action Tr[i 7 ' 2 ly ]. 
There is no spin current from the non-Abelian gauge field 
when it is assumed to be the frozen background field.) 

When the magnetic ordering occurs, the wavefunctions 
are further projected onto the spin component at each 
site, which is described by the field operator ip a of elec- 
trons decomposed into 



?/v(x) = z CT (x)/(x) 



(14) 



with / being the spinless fermion corresponding to 
the charge degrees of freedom while z a being the 2- 
component spinor corresponding to the direction of the 
magnetization at x. Putting ea. (fl4l into eq. p^)) . we ob- 
tain the effective Lagrangian for /-field as 



id - 



so 



■A + 



(V 



ieA) 2 



2m 



where 



i(z\d fl \z) is the U(l) field originating from 
the Berry connection of the spin wavefunctions, and 
af° = A a a {z\T a d a \z), = Ag(z|r a |z) are U(l) field 
originating from the SOI. 

Here a remark about t he g eometrical meaning of the 



gauge field is in order 14 1. The magnetic flux made 



from a B , i.e., b = V x a B , corresponds to the solid angle 
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subtended by the spins. Namely, the scalar spin chiral- 
ity given by Si • (Sj x S fe ) is reduced to b along the 
direction normal to the plane made by the 3 sites i, j, 
k in the continuum limit. Therefore, the effective mag- 
netic field is produced when the spin structure is non- 
coplanar. In summary, there are three "electromagnetic 
fields" in magnetic systems, i.e., (i) a B due to the Berry 
phase associated with the non-coplanar spin structure, 
(ii) a so from the spin-orbit interaction (SOI), and (hi) 
usual Maxwell EMF A. 



C. Emergent electromagnetic field in correlated 
electronic systems 

Up to now, we have considered the gauge fields in the 
band structure or magnetically ordered state, i.e., the 
single-particle properties, where the gauge fields arc the 
static background fields. The electron correlation effect 
corresponds to the fluctuating spin beyond the mean field 
theory, and accordingly the gauge fields become dynami- 
cal [151 ] . This issue has been studied extensivelyin the re- 
search of high temperature superconductors [l4|, IH| • As 
the simplest example, let us consider the non-linear sigma 
model describing the low energy physics of the quantum 
antiferromagnet 0, [I?} • 

S=-J^ dr J dr\d^n\ 2 , (16) 

where the unit vector n specifies the staggered moment, 
and g oc 1/5 (S : spin quantum number) is the dimen- 
sionless coupling constant representing the quantum fluc- 
tuation. We introduce the spinor field z =* (z^,z±) with 
the constraint Ylcr z t z a — 1> n can be represented by 
n = z\(o) a a i z a i with a = (a x , a y , a z ) are the Pauli ma- 
trices. Then S in eq. (|16[) can be written as 

S=-\(d fl +ia fl )z <7 \ 2 (17) 
9 

with the gauge field = —iz\,d^z a corresponds to the 
Berry phase connection between the 2 spinors at the 
neighboring 2 spatial points. Therefore, this corre- 
sponds to the U(l) gauge field due to the spin chirality 
discussed in the previous subsection. When we regard 
Zfj as being the classical field, i.e., the condensed compo- 
nent corresponding to the magnetic ordering, the gauge 
field a M is reduced to a B discussed in the last subsec- 
tion. Without the magnetic ordering, on the other hand, 
there is no condensation of z-field and hence the Higgs 
phenomenon for the gauge field does not occur. In 
this spin liquid state, this gauge field is the dynamical 
degrees of freedom, i.e., obtains the effective action by 
integrating over the rapid components of the spin fluc- 
tuations (Note that there is no term describing the 
dynamics of the gauge field at the starting since it is in- 
troduced to express the constraint.) Therefore the gauge 
field appears dynamically as emergent phenomenon. The 



quantum spin fluctuations and electron correlation ef- 
fects are the subjects of intensive studies in condensed- 
matter physics, and the gaug e theoretical approach has 
been a powerful tool [1, IB1 Il6l| . Historically, motivated by 
the studies on strongly correlated systems and quantum 
Hall systems from the gauge theoretical viewpoint dur- 
ing 1990's, the conventional materials, which can be well 
described by band theory and/or the mean field theory, 
have been revisited from this new perspective. Below, 
we take some of the examples, where new physical ef- 
fects have been explored in conventional systems from 
the viewpoint of gauge fields. 



II. MATERIALS AND PHENOMENA OF 
EMERGENT ELECTROMAGNETISM 

As described above, there are several sources of the 
EEMF, and correspondingly there are so many related 
phenomena in condensed-matters. The common feature 
of these phenomena is that the topological currents play 
an essential role. In solids, there are many imperfec- 
tions such as the defects and impurities, which cause the 
scattering of electrons and dissipation. This dissipation, 
i.e., Joule heating, is balanced with the energy supplied 
from the external electric field. This is the usual Ohmic 
current. In the superconducting state, the macroscopic 
quantum coherence and the associated "rigidity" against 
the phase twist by the external magnetic field produces 
the superconducting current. This current does not de- 
cay nor cause any dissipation since the current-flowing 
state is in the thermal equilibrium, and is separated by 
the macroscopic energy barrier from the zero current 
state. To this list we would add here the third cate- 
gory of the current, i.e., the topological current, induced 
by the gauge field or curvature of the Hilbert space, i.e., 
EEMF. This third category includes the quantum Hall 
current and polarization current in ferroelectric materi- 
als, which are related to the Berry phase. This topolog- 
ical current does not require the off-diagonal long range 
order (ODLRO), and can exist in the normal states even 
at room temperature. Therefore, the topological currents 
will be of vital importance when the applications to elec- 
tronics are considered. We will describe below several 
novel phenomena driven by this topological current in 
solids. 



A. Multiferroics 

The minimal coupling to the SU(2) gauge field A® in 
eci.([T2"j) or in eq. (fT5)) for magnets means the close relation 
between the spin current and the electric polarization 
P. More explicitly, P is given by the derivative of the 
Lagrangian with respect to E, i.e., 

P l oc e iav fi, (18) 
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which means that the spin current produces the ferro- 
electric moment. Then the next problem is how one can 
produce the spin current. In the magnetically ordered 
state, the expectation value of the spin current 0„)) in 
the ground state can be nonzero for the noncollinear spin 
configuration. Consider the expression of the spin cur- 
rent operator in the tight-binding model 



3% = i^2t rl +^c\ a ((J a ) a0 c l+f _ l p + h.c, 



(19) 



where tij is the transfer integral, a a (a = x, y, z) the Pauli 
matrices, and a, /3 the spin indices. For simplicity, let us 
consider the case of the spin ordering within the iry-planc, 
and the corresponding wavefunciton Xi =t ("Tf'l/f")- 
Taking the expectation value with this wavefunction as- 
suming the half-filling, one obtains 



°ii-\-fi 



sin 



(4>i - i+M ), 



(20) 



where we employed the perturbation theory in t/U (U: 
on-site Coulomb interaction or the Hund's coupling). 
This is analogous to the Josephson current in super- 
conductor with (j> being the phase of the order param- 
eter, and the spin supercurrent in magnet is induced 
by the tilted spins. Putting this into ea. ([T5| . (P z ) cx 
£znv sm (0i — 4>i+v)- More general expression is easily ob- 
tained for the polarization P^ obtained by the 2 spins S, 



and Sj as 



P = Tjetj x (Si x Sj) 



(21) 



with r\ being the coupling constant related to SOI, and 
e,;,- is the unit vector connecting the 2 sites i and j [HI, 
19j . The vector product of the 2 spins Xij = Sj x Sj is 
called vector spin chirality. This quantity is even with 
respect to the time-reversal operation T similarly to the 
spin current. 

This generic argument is applied to the real materials 
as follows. Compared with the free electrons in vacuum, 
the strength of the relativistic spin-orbit interaction can 
be enhanced by the factor of ~ 10 6 which is the ratio 
of the rest mass of the electrons mc 2 to the band gap. 
For 3d electrons in transition metal atoms, the SOI A is 
typically of the order of ~ 20 — 50meV, while it becomes 
~ 0.5eV for 5d electrons. The electron correlation en- 
ergy, on the other hand, decreases from 3d to 5d since 
the wavefunction is more and more expanded. In the 
cubic crystal field in transition metal oxides, the 5-fold 
degeneracy of d-orbitals is lifted due to the ligand field 
of oxygens. As the result, 3-fold degenerate i 2ff orbitals 
(xy, yz, zx-orbitals) with lower energy, and doubly degen- 
erate e g -orbitals (x 2 — y 2 ,3z 2 — r 2 -orbitals) with higher 
energy are formed. The matrix elements of the orbital an- 
gular momentum £ are zero within the e g -orbitals. On the 
other hand, they are nonzero among i^g-orbitals and also 
between the e g - and i 2 g-orbitals. This SOI is the origin 
of the relativistic coupling between the magnetism and 



electric polarization. To give a more explicit prediction 
for transition metal oxides, the cluster model of mag- 
netic ions sandwiching an oxygen ion has been studied 
theoretically by taking into account the spin-orbit inter- 
action when deriving the super-exchange interaction [l8| . 
As mentioned above, the spin current flows between the 2 
non-collinear spins Si and Sj , which produces the electric 
polarization P as given by 

p " ~¥ ( a" } Ieij x (Sl x Sj) ' (22) 

where / = (p x \z\d zx ), and A (V) is the energy difference 
(hybridization) between the p-orbitals and the d-orbitals. 
The SOI interaction is implicitly included in this model 
by picking up one doublet after splitting by the SOI. Ap- 
plying this result to various magnetic structures, one can 
easily predict the presence or absence, and the direction 
of the polarization. This theory is not contradicting with 
the symmetry argument developed for magnets [20I [2l| , 
but stresses the physical mechanism of the spin current- 
induced polarization. 

A recent experimental breakthrough is the discovery of 
the multiferroic behavior in PMn0 3 (P=Cd,Tb,Dy) @, 
3]. In these materials, the spontaneous electric polariza- 
tion P s appears accompanied by the magnetic order, and 
they are necessarily strongly coupled. It has been also 
revealed that the magnetic structure which induces the 
electric polarization is the cycloidal spiral in good ac- 
cordance with the theoretical prediction above [22T - [2~ij ] . 
Figure 4 shows the representative multiferroic behavior 
of DyMn03. In this material, the spin rotation plane flop 
from be- to ab-planes, and accordingly the direction of 
the electric polarization changes from c- to a-axes [![. 

Now the extensive experimental studies have been done 
to explore the multiferroic materials, ejj., NI3V2O8 251. 
Ba .5Sri^Zn 2 Fei 2 O22 SI, CoCr 2 4 [13, MnW0 4 H, 
CuFe0 2 ^i, LiCuV0 4 [3^, andLiCu 2 2 ^ are discov- 
ered to be multiferroics. Namely, multiferroicity is not 
a special phenomenon but is a rather ubiquitous phe- 
nomenon in Mott insulators. These experimental find- 
ings urged the systematic theoretical studies on the mi- 
croscopic mechanisms of the spin-related electric polar- 
ization [H, [33| . The perturbative approach in both V/ A 
and A/ A is employed, where V and A represent the trans- 
fer integral and the charge transfer energy between the 
transition-metal (TM) d and ligand (L) p orbitals. This 
analysis concludes that the polarization P r +s appearing 
on the bond between the sites r and r + e is given by 

P r+ | = P ms (m r • m r+e )e + P sp e x (m r x m r+e ) 

+ P orb [(e • m r )m r - (e • m r+e )m r+e ] , (23) 

where m r is the spin direction at r. The first term 
pms ^ (y/A) 3 is the polarization due to the exchange- 
striction, which is nonzero when the inversion symmetry 
between r and r + e is absent because the 2 intermediate 
states of doubly occupied d-orbitals becomes inequiva- 
lent. This term does not requires the SOI, hence is con- 
sidered to be larger than the rest of the terms if it exits. 
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The second term P sp oc (A/A)(V/A) 3 is due to the spin 
current mechanism already discussed. The third term 
P orb ~ min(A/V, 1)(V/A) is nonzero for the partially 
filled t2 g orbitals and comes from the modification of the 
single-spin anisotropy due to the electric field [H, [33j . 
These three contributions appear differently depending 
on the wavevector of the polarization. Now the origins of 
the multiferroic behavior of various materials are investi- 
gated and are classified according to these 3 mechanisms 
in most of the cases 0. 

In eq. (|23]l . the polarization is due to the relativistic 
SOI (spin current mechanism or spin anisotropy mecha- 
nism) or the broken inversion symmetry at the center of 
the bond (exchange-striction). Therefore, it is an impor- 
tant question what is the most ubiquitous and generic 
interaction between the electric field and the spin system 
without these mechanisms. By a careful analysis of the 
super-exchange processes in Mott insulators taking into 
account the quantum dynamics of the spins in the in- 
termediate states, we have derived the following effective 
Lagrangian [34| 

L E = - J d 3 rT ab E a e b , (24) 

where E is the electric field and 

e = i sin 9(d t eVtp - d t ipV6) (25) 

is the electric field component of EEMF associated with 
the spin structure derived from the U(l) Berry connec- 
tion. Here n = (cos (p sin 6, sin tp sin 9, cos 9) is the unit 
vector of the magnetization in the continuum approxima- 
tion, and T a b is the tensor which depends on the details of 
the atomic configurations. One important remark here 
is that T a b vanishes in the single-orbital case since the 
spins form singlet and has no quantum dynamics in the 
intermediate states in the processes of the super-exchange 
interaction. However, in most of the cases, the orbital de- 
grees of freedom is active and T ab is nonzero. This cou- 
pling is analogous to the spin motive force in the metallic 
ferromagnetic systems by which the domain wall or vor- 
tex motions produce the voltage drop [35U37j | . 



B. Topological Hall effects 

As explained in eq.©, the most natural phenomena 
expected from the gauge fields in momentum space are 
related to the transverse motion of the electrons to the 
external electric field, i.e., the Hall effects. Contrary to 
the usual Hall effect driven by the Lorentz force, we call 
the Hall effects originating from the gauge fields as topo- 
logical Hall effects, and we describe below some of the 
examples. 



1. Spin-orbit interaction and anomalous Hall effect 

Hall discovered the Hall effect in metallic ferromag- 
netic systems due to the spontaneous magnetization in- 
stead of the external magnetic field (38|. This effect is 
called anomalous Hall effect (AHE) and its mechanism 
has been a controversial issue for more than a century. 
It is agreed that the effect is due to the SOI combined 
with the spontaneous magnetization, but the issue was 
the role of the impurity scatterings. Intrinsic mechanism 
was first proposed by Karplus-Luttinger (KL) [39(, who 
considered the anomalous velocity for the first time dis- 
cussed in section I. However, their theory has been crit- 
icized by Smit |4(| saying that the impurity scattering, 
which is inevitable and also indispensable to reach the 
steady state under electric field, invalidates the KL the- 
ory. This posed the question if the dissipationless topo- 
logical current can survive even in the presence of the 
dissipativc Ohmic current. Smit [40jj proposed the skew 
scattering mechanism, where the impurity scattering in 
the presence of SOI gives the asymmetry of the transi- 
tion rates between Wk->-k' and Wk'->k- This breakdown 
of the detailed balance leads to the net current perpen- 
dicular to the applied electric field, i.e., the Ohmic trans- 
port current is slightly distorted in the perpendicular di- 
rection. Later Berger (4f| proposed another mechanism 
called side-jump, where the transverse shift of the elec- 
tron trajectory occurs at the scattering in the presence of 
the SOI. These 2 mechanisms are called extrinsic ones as 
opposed to the intrinsic one by KL. The recent advances 
in this problem is two-folds [38]. One is the recognition 
that the band structures of the ferromagnetic materials 
can be topologically nontrivial characterized by Chern 
numbers. Haldane [42j was the first to notice that the 
quantum Hall effect can be realized without the Lan- 
dau level formation in a tight-binding model on a lattice 
with the complex transfer integrals. What is recognized 
recently is that this scenario can be realized in the ferro- 
magnetic materials with the SOI [43]. In this case, even 
though the overlaps of the band dispersions make the sys- 
tem metallic, one can define the Chern number of each 
band when the gap opens at each k-point. Therefore, the 
ferromagnetic metal could be an implicit quantum Hall 
system with A: z -dependent Chern numbers. The other is 
the state-of-art first-principles calculations of electronic 
states and the Hall response as well as the accurate ex- 
perimental measurements of the AHE in various materi- 
als [38|. In short, the transverse conductivity a xy due to 
the intrinsic mechanism can be written in terms of the 
Berry phase curvature, i.e., the magnetic flux 6 n2 (k) for 
the Bloch wave state with crystal momentum k and the 
band index n [Hj]: 

nk 

where the /(e n (k)) is the Fermi distribution function for 
the energy e„(k) of the Bloch state, and the vector poten- 



7 



tial a njU (k) is defined in eq.([2]), and b„(k) = x a„(k). 
The meaning of this expression is that the Hall current is 
the sum of the anomalous velocities of Bloch states occu- 
pied by electrons in the equilibrium distribution. There- 
fore the transverse conductivity a xy represents the gauge 
field distribution in the momentum space. Especially 
when the chemical potential is in the gap, the k-integral 
in eq. (|2l)| is over the whole first Brillouin zone (1st BZ), 
and one might think it is zero due to the periodicity with 
respect to k — > k + G with G being the reciprocal lattice 
vector. Namely the contour integral over boundary of 
the 1st BZ appears to cancel. However this is not always 
the case, and the integral is (e 2 /h) x integer especially 
in the 2-dimensions. This integer is the topological num- 
ber called Chern number for the U(l) fiber bundle of the 
Bloch wavefunction (44j . The finite Chen number means 
that the phase of the Bloch wavefunction can not be de- 
fined continuously over the whole 1st BZ analogous to 
the Yang-Wu construction of Dirac monopole |4j, [5j, |45j . 
Here the plural patches need to be introduced to cover 
the 1st BZ. The relevance of the Dirac monopole to the 
Hall effect is understood more explicitly in the following 
3-dimensional model. 

H(k)= Y, k ^a- (27) 

The corresponding magnetic field associated with the 
Berry phase reads 

b ± (k) = ±^3, (28) 

where ± corresponds to the 2 bands with energy e±(k) = 
±|k|. This means that the magnetic monopole exists at 
the band crossing point k = 0, and Vfc-6±(k) = ±2-7T(5(k) 
has the delta-functional singularity. The contribution 
from the fixed k z is given by 

e 2 

cr xy (k z ) = ±— sign/c 2 . (29) 

For the 2-dimensional system where k z is regarded as a 
parameter m characterizing the system, it describes the 
quantum Hall system with a xy = ±|^. The factor 1/2 
is not realized in real system because the Dirac fermion 
appears always in the pair, i.e., species doubling, in the 
1st BZ and hence the Chern number is an integer. 

Essential to AHE is the fact that the nontrivial topo- 
logical structure described above is not special as in the 
case of quantum Hall system, but is generic for mag- 
nets in the presence of the SOI. A simple 2-dimensional 
3-band model for ti g orbitals was constructed, and the 
uniform magnetization produces the nonzero Chern num- 
ber for each band |43j. In the absence of the SOI, the 
up and down spin bands are independent of each other 
except the exchange energy splitting, and there occurs 
several band crossing points but no finite Chern num- 
bers. Then SOI lifts the degeneracies to generate the 
mass term m as discussed above and produce the finite 



Chern numbers. This is the reason why the SOI can not 
be treated perturbatively in sharp contrast to the previ- 
ous theoretical treatments. This nontrivial behavior has 
been found both in the first-principles band calculation 
and in the experiment in SrRu03, a metallic ferromag- 
net with T c = 130K [46]. In Figs. 5 shown the transport 
properties of SrRu03. Figure 5(c) shows the temperature 
dependence of the Hall resistivity which is dominated by 
the anomalous contribution. Just below T c , it is negative 
and turns into positive and again has the maximum. This 
characteristic behavior strongly suggests that the pertur- 
bative expansion in AM (A: SOI, M: magnetization) is 
not allowed and the AHE is a fingerprint of the Berry 
phase of the Bloch wavefunction. We plot in Fig. 5(d) the 
transverse conductivity a xy as a function of the sponta- 
neous magnetization compared with the first-principles 
band structure calculation. This non-monotonous tem- 
perature dependence including the sign change is due to 
the band crossings acting as the magnetic monopoles in 
momentum space. Accordingly Fig. 6 shows the dis- 
tribution of the Berry curvature b nz (k) as a function of 
(k x ,k y ) at fixed k z — 0. It shows a sharp peak around 
k x = k y — and the ridges along k x = ±k y . This sharp 
peak represents the monopole corresponding to the band 
crossing. When the Fermi energy ef is near the magnetic 
monopoles, the electrons are subject to the strong gauge 
field and hence the contribution to the Hall constant is 
resonantly enhanced. The integral over k z leads to the 
partial cancellation of the positive and negative contri- 
butions, but still the rapid change of a xy as a function 
of Ef results. This explains both the first-principle band 
calculation and experimental results, which are shown in 
Figs. 5(c) and (d). Similar conclusion is obtained also for 
the AHE in Fe [47J, where the sharp spiky structures of 
the Berry phase distribution has been found. 

An important recent development is the experimental 
measurements of dynamical AHE, i.e., <j xy (uj) in ferro- 
magnetic materials. Thanks to the rapid progress of the 
telahertz spectroscopy [H, |49| , the low frequency struc- 
tures of the dynamical AHE are being revealed. Espe- 
cially, the dominance of the contributions from the band 
crossings is expected to appear in the w-dependence of 
a xy (ui) as well as the non- monotonous temperature de- 
pendence of a xy {ui — 0). Figure 7 shows an example 
of the low frequency data for AHE in SrRu03. It is 
clearly seen that nontrivial structure exists at around 
4meV, which is well fitted by a model assuming the band 
crossing near the Fermi energy. 

Up to now, it is argued that in some materials, the 
AHE can be explained by the intrinsic mechanism based 
on the Berry phase concept, but the controversy between 
the intrinsic and extrinsic needs to be resolved. As for 
this problem, again the role of the enhanced Berry cur- 
vature near the band crossings and the topological na- 
ture are crucial. In ref. 50] a model with a band anti- 
crossing is considered and its Hall response is calculated 
taking into account the disorder scatterings. The Hall 
conductivity a xy is obtained as a function of a xx repre- 
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senting the strength of the disorder. Figure 8 summa- 
rizes the results including the various experimental data 
plotted. (Note that the calculation has been done for 
a 2D model, and cr's are scaled by the lattice constant 
along the z-direction which is assumed to be = 0.4nm.) 
This result captures the respective role of the transport 
and topological currents. In the very clean region where 
a xx > 10 ft cm , the skew scattering mechanism (ex- 
trinsic mechanism) is dominant and hence a xy oc cr xx . 
When the disorder strength becomes larger, the trans- 
port current is suppressed much more rapidly compared 
with the topological current, the latter of which is "pro- 
tected" topologically, leading to the plateau region in 



the region 10 ft 



n 1 cm 1 < Oxx 



< io 6 ^-^ 



In this 



intermediate region, to which many of the metallic fer- 
romagnets belong, the intrinsic mechanism is dominant 
and hence the first-principles calculations can predict the 
AHE semiquantitatively. When the disorder becomes 
even stronger, a xy begins to decrease obeying the ap- 
proximate scaling law a xy oc cr^ 6 . This scaling is in good 
accordance with the experimental results as seen in Fig. 
8, but its understanding is still left for future studies. To 
conclude this subsection, AHE offers an interesting lab- 
oratory where both the Ohmic transport and topological 
currents coexist, and it has been established that the lat- 
ter can contribute to the Hall effect even in the metallic 
systems. Quantized version of the AHE, where most of 
the bulk states are localized [5l[ or the chemical poten- 
tial is in the gap (52j, have been discussed also. In this 
case, the current is carried by the edge channels as in the 
case of quantum Hall system. 



2. Scalar spin chirality and anomalous Hall effect 

There are a vast variety of the non-collinear spin struc- 
tures found experimentally. These systems are the ideal 
laboratory to test the idea of U(l) EEMF discussed in 
section I. Of particular interest is the idea of "scalar spin 
chirality" Xijk defined for the 3 spins as 



Xijk — Sj • [Sj X S/j), 



(30) 



which is T-odd. As discussed in section IB and Fig. 2, 
the scalar spin chirality acts as an effective magnetic field 
and leads to the Hall response [53,, 54]. First, we consider 
the case where the magnetic unit cell is small, and the 
Bloch wavefunctions are defined in the 1st BZ. The total 
gauge flux penetrating the unit cell is zero or integer 
multiple of 2ir due to the periodicity since the contour 
integral of a(r) along its boundary vanishes ( mod 2tt) 
due to the periodicity. On the other hand, the gauge 
field corresponding to the Berry phase becomes nonzero 
for the Bloch wavefunctions in the momentum space with 
a possible finite Chern number of each band. Roughly 
speaking, this happens when the unit cell contains the 
more than 2 different types of the loop, and each band 
feels the fluxes with different weight to obtain the Chern 
number. Ohgushi et al. [54| showed that this idea can 



be realized in a ferromagnet on Kagome lattice with the 
non-coplanar spin configuration. 

Pyrochlore lattice can be regarded as the 3-dimensional 
generalization of the Kagome lattice, i.e., pyrochlore lat- 
tice contains the Kagome lattice normal to the [111] di- 
rection. In the material Nd2Mx>207 (hereafter we de- 
note it as NMO), there are 2 interpenetrating sublat- 
tices composed of the tetrahedrons of Nd and Mo atoms 
shifted along the c-axis as shown in Fig. 9(a) 55]. The 
conduction electrons are on the Mo sublattice, while the 
localized spins on Nd sublattice are subject to the strong 
anisotropy to form non-coplanar spin configurations. The 
localized spins and conduction electron spins are coupled 
antiferromagnetically, which transmit the scalar spin chi- 
rality of Nd to conduction electrons (Fig. 9(b)). Figure 
9(c) shows the temperature dependence of the Hall resis- 
tivity which shows the steep increase as the temperature 
is lowered corresponding to the increase of the exchange 
field from Nd spins and associated scalar spin chirality. 
The absolute value of the low temperature value is consis- 
tent with the theoretical calculation taking into account 
the tilting angle of the Mo spins estimated from the neu- 
tron scattering experiment [551 ] . Furthermore, the obser- 
vation of the sign change of a xy under the external mag- 
netic field along [111] direction [5(| is consistent with the 
sign change of the spin chirality. 

The analysis in k-space above is justified when the 
magnetic unit cell is small compared with the mean free 
path I. In this case, the EEMF is defined by the Berry 
phase in k-space. In the other limit, i.e., when I is shorter 
than the size of the spin texture £, the EEMF is defined 
more appropriately in r-space, which acts as the effective 
EMF with uj c t << 1 (uj c : cyclotron frequency, t: mean 
free time). This situation is realized in the recently dis- 
covered Skyrmion crystals in noncentrosymmetric mag- 
nets with B20 structure [13, El]. MnSi, (Fe,Co)Si, and 
MnGe are the examples of this class of materials, which 
has Dzyaloshinskii-Moriya (DM) interaction as described 
by the following Hamiltonian in the continuum approxi- 
mation. 



H = dr 



^(VM(r)) 2 +aM(r)-(VxM(r))-H-M(r) 

(31) 

where M(r) is the magnetization, H the magnetic field, J 
the ferromagnetic interaction, and a the DM interaction 
constant. The lattice constant is taken to be unity. The 
ground state under zero magnetic field is the single spiral 
with the spins rotating in the plane perpendicular to the 
wavevector q. The length of q = |q| is determined to 
the the ratio a/ J while its direction is determined by the 
weak spin anisotropies not included in the Hamiltonian 
eq.dH]). 

Recently, a neutron scattering experiment identified 
the mysterious A-phase in MnSi as the Skyrmion crystal 
state stabilized b y th e external magnetic field and ther- 
mal fluctuations [5 71 ]. The conical spin structure is the 
most stable state under the magnetic field in the ma- 
jor part of the phase diagram in 3D crystal. On the 
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other hand, when one reduces the thickness of the sam- 
ple smaller than the wavelength of the spiral, the coni- 
cal state becomes energetically higher than the Skyrmion 
crystal when the external magnetic field is perpendicu- 
lar to the film. Actually, a Monte Carlo simulation of 
the 2D magnet with DM interaction concluded that the 
Skyrmion crystal state is stable in much wider region of 
the phase diagram in (T, B)-plane (T: temperature. B: 
magnetic field ) including the zero temperature case[59|. 
Motivated by this study, a recent experiment using the 
Lorentz tunneling electron microscopy (TEM) succeeded 
in real-space observation of the Skyrmion crystal in the 
thin film of (Fe,Co)Si [HI]. Figure 10 shows the phase 
diagram and the Lorentz TEM images of the Skyrmions 
and Skyrmion crystals. The left panel shows the experi- 
mental results, while the right panel the results by Monte 
Carlo simulations, which indicate the good agreement. 

As shown in Fig. fl schematically, each Skyrmion 
wrap the sphere once in the order parameter space of M, 
which means that the integral of the magnetic flux b || z 
over one Skyrmion is 2n. From this one can estimate 
the effective magnetic field induced by the Skyrmion. 
It is typically 4000T when the size of the Skyrmion 
£ = Inm, and inversely proportional to the square of 
£ 2 . Since MnSi, (Fe,Co)Si, MnGe, are metallic systems, 
the conduction electrons are coupled to the spins and 
spin chirality. Therefore, the Hall effect is expected due 
to the Lorentz force driven by b replacing the exter- 
nal magnetic field B. Recently, the Hall measurement 
on MnGe has been done and the contribution from the 
Skyrmion has been analyzed [6(|. The Hall resistivity 
Ap yx = 200nf2cm is due to this topological contribu- 
tion, which is corresponding to the effective magnetic 
field b = HOOT for £ = 3nm. These values are com- 
pared with the case of MnSi, where £ = 18nm, b = 28T, 
and A.p yx = 5nf2cm. 

The dynamics of Skyrmions and Slyrmion crystal is an 
intriguing issue to be studied. A recent experiment [6l| 
shows the Skyrmion crystal motion driven by a current 
with the threshold value j c = 10 2 A/cm 2 m, which is or- 
ders of magnitude smaller than that for the domain wall 
motion in ferromagnets [HJ . Theoretically, the motion of 
the Skyrmions will induce the electric field e as electro- 
motive force due to dh/dt. This leads to several physical 
consequences such as a new mechanism for the damping 
of the Skyrmion motion, and the motion of the Skyrmion 
transverse to the current (Skyrmion Hall effect) |63f . 



3. Hall effect of light 

The topological Hall effect does not require that the 
particles are charged, and hence is possible even for the 
uncharged particles. A representative example of un- 
charged particle is the photon, and one can ask if the 
Hall effect of light can occur or not. The constraint which 
induces the Berry curvature in this case is that the elec- 
tric and magnetic fields of light is always perpendicular to 



the direction of the propagation, i.e., that of the wavevec- 
tor k. This situation corresponds to the strong coupling 
limit of SOI in electrons. Actually, the effect of the Berry 
phase in the propagation of light in the optical fiber has 
been studied both experimentally and theoretically, and 
the rotation of the polarization was confirmed for the 
light with the k- vector slowly changing and enclosing an 
area on the sphere (64| . 

Note that the quantum nature of the photon is not re- 
quired for this Berry phase effect, which is solely due to 
the wave nature of light. Therefore, this offers a unique 
opportunity to study the Berry phase effect in classical 
systems. Extending this idea, we have derived the semi- 
classical equation of motion for the wavepacket of light 
including the effect of the finite wavelength (65[. (Note 
that the wave-optics is reduced to the geometric optics 
in the short- wavelength limit, the former of which corre- 
sponds to the quantum mechanics while the latter to the 
classical Newtonian mechanics.) This semiclassical equa- 
tions of motion for the position r, the wavevector (mo- 
mentum) k and the spin state \z) (2-component spinor 
describing the polarization of light) reads 



dt 
dk 
~dl 

dt 



— = -[V«(r)]k, 



(32) 



where v(r) — c/n(r) is the velocity of light in the 

medium, [AJaa' = — ie^ k VkeA'k is the connection of 
2x2 matrix form, i.e., SU(2) gauge connection, and 
J^k = Vk x Ak + Ak x Ak is the field strength. It is 
diagonal in the basis of the right and left-circular polar- 
ization, i.e., Jlk = C3k/|/c| 3 . This equation indicates that 
the Berry curvature has the opposite sign for opposite 
circularly polarized light, which gives the polarization- 
dependent anomalous velocity adding to the group ve- 
locity of wavepacket. More explicitly, one can consider 
the reflection and transmission of light at the interface 
of 2 media with different dielectric constants [65]. The 
anomalous velocity is finite when the particle is subject 
to the gradient of the refractive index Vn(r), and results 
in the finite transverse shifts of the transmitted and re- 
flected lights at this reflection/transmission. These shits 
are in the opposite direction for opposite circular polar- 
ization. This effect has recently been observed experi- 
mentally using the "weak measurement" 66]. 

The idea of the Berry phase for light can be general- 
ized to the photonic crystal, where the Bloch waves are 
formed. Especially, one can consider the crystal with 
distortion, where the Berry phase is defined in the 6- 
dimensional space (k, r). Especially, the curvature is en- 
hanced near the gap edge when the gap is small. This 
situation is realized for the X-ray propagating in crystals 
since the deviation of the dielectric function from unity 
is typically of the order of 10 -6 , and hence the "periodic 
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potential" for the Bloch wave is very weak. Sawada et 
al. [§3] considered the role of this enhanced Berry cur- 
vature in the propagation of X-ray in deformed crystal. 
The prediction is that the shift in the trajectory Ar of 
the X-ray is given by 



(33) 



where G is the reciprocal lattice vector satisfying the 
Bragg condition for the wavevector k of the X-ray, u 
the displacement of the crystal, to the frequency of X- 
ray, and Aw is the gap in the Bloch wave of X-ray. As 
mentioned above, lj/Auj can be as large as 10 6 , which de- 
termines the magnification of Ar compared with u con- 
sidering that |G| ~ |k|. This theoretical prediction has 
been recently confirmed experimentally using the single 
crystal of Si [68| , and offers a new principle for the X-ray 
microscope. 



4- Hall effect of magnons 

Another example of the uncharged particle is the 
magnons in insulating magnets. Although the charge 
current is not available in this case, the energy current 
can be carried by the magnons and hence the thermal 
Hall effect without resorting to the Lorentz force is When 
one considers the propagation of the magnon wave in a 
ferromagnet, it behaves as a Bloch wave and the band 
structure is formed. When more than 2 atoms are in the 
unit cell, the Berry phase is generally expected for this 
Bloch wavefunction. We consider the Hamiltonian with 
the Dzyaloshinskii-Moriya (DM) SOI as 

h = -J2 J ^ ■ s j + B i ■ ( s « x s j) - £ H ■ s - ( 34 ) 



Suppose the external magnetic field H and ferromagnetic 
moments align along +z-direction. Then, the spin wave 
Hamiltonian can be written as 

H sw = ~Y, 4-#(e~^ w &l&j + h.c.) + H z b\b u (35) 



where Jije il *' ij = Jij 



iD?,. This indicates that the DM 



interaction acts as a vector potential and effective mag- 
netic flux for the propagating magnons [69j . Similar to 
the case of AHE, the distribution is such that the to- 
tal effective magnetic flux is zero (or equivalently the 
integer multiple of 2tt) when integrated over the unit 
cell. Therefore, finite effect survives when the inequiv- 
alent loops exist in the unit cell. Kagome lattice and its 
3D generalization pyrochlore lattice are the representa- 
tive crystal structures satisfying this condition. And ac- 
tually the thermal Hall effect has been recently observed 
in an insulating ferromagnet LU2V2O7 with pyrochlore 
structure |69| . In the pyrochlore structure, the midpoint 
between any 2 corners of a tetrahedron is not an inversion 



symmetry center, and hence the nonzero DM interaction 
is expected. Symmetry further determines the direction 
of the DM vectors Dy as shown in Fig. 12(a). 

In this material, the spontaneous magnetization M 
emerges below Curie temperature T c — 70K, which is 
isotropic and almost coincides with 1 Bohr magneton 
(Pb) at low temperatures, indicating the collinear fer- 
romagnetic state with spin S = 1/2. Below T c , the 
thermal Hall conductivity k,r is discernible, whereas it 
is very small above 80 K. The magnitude of the ther- 
mal Hall conductivity has a maximum at around 50A'. 
Similar to the magnetization, the thermal Hall conduc- 
tivity steeply increases and saturates in the low magnetic 
field region. Therefore, the observed thermal Hall effect 
is due to the spontaneous magnetization M as in the 
case of AHE. kh gradually decreases with increasing the 
magnetic field furthermore due to the opening of the gap 
in the magnon spectrum induced by the magnetic field. 
All those effects are contained in following theoretical ex- 
pression for Ka/3 [6!| : 



pi 



3 / 2 ha 



k B T 
2JS' 



Li 



5/2 



exp 



9PbH \ 
2JS J 

k B T J 



(36) 



where <& a p — — e Q( 3 7 n 7 _D/(8v / 2 J) ( £ Q( g 7 : totally an- 
tisymmetric tensor, n: direction of the magnetiza- 
tion, a: lattice constant, D = |Dy|), and Li„(z) = 
Sfe=i 00 z Jk n . Figure 12(b) shows the comparison be- 
tween the experimental result and ea. (!36p . The only ad- 
justable parameter is the ratio D/J = 0.32, which is a 
reasonable value for transition metal oxides. From these 
facts, it is convincing that the thermal Hall effect ob- 
served in LU2V2O7 is due to the magnons affected by the 
Berry curvature due to the DM interaction. 



5. Spin Hall effect 

The concept of the anomalous velocity can be gener- 
alized to the time-reversal T— symmetric systems. As 
already discussed in eg. (|3"2"1) . the SU(2) Berry curvature 
is non-zero even in the T— invariant system, and the 
analogous effect is expected for the electrons also. This 
idea has been explored for the band structure of semi- 
conductors [7(j. In GaAs and Ge, the valence bands 
consist of three p-ortbitals (L = 1) while the conduc- 
tion bands s-orbitals (L = 0). Therefore the 6 bands 
(3 times the spin degeneracy 2) constitute the valence 
bands, which are split into 2-fold degenerate split-off 
(SO) bands (J = 1/2), and the 4- fold degenerate bands 
(J = 3/2) at the T-point (k = 0). The 4-fold degen- 
erate bands are further split into 2 doubly-degenerate 
heavy-hole and light-hole bands, i.e., 2 Kramers degener- 
ate bands at finite k. (Here we neglect the breaking of the 
inversion symmetry in the case of GaAs, which is small.) 
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These 4 bands are described by Luttinger Hamiltonian 

as 



^(k) = ^ ((7i + ^72)fc 2 - 2 72 (k • S) 2 j .(37) 

As discussed above, the degeneracy acts as the mag- 
netic monopole for the Berry curvature, and the T-point 
is especially interesting since the 2 doubly-degenerate 
bands touch there. Since there are Kramers degeneracy 
at each k-point, the Berry connection is the 2x2 matrix 
(SU(2)). There are 2 pseudo-spin states, i.e., parallel and 
anti-parallel pseudo-spin to k = k/|k|, which is called he- 
licity. The gauge held and hence the anomalous velocity 
depend on the helicity of the state, i.e., it has the op- 
posite sign for different helicity states. Therefore, even 
though there is no net charge current, the external elec- 
tric field can induce the transverse spin current. This is 
called the (intrinsic) spin Hall effect, and for the case of 
the p-doped GaAs/Ge the spin current jg was predicted 
to be 



dop< 



J) = a s H e» k E k 



with 



' H 



6tt 2 



{hp — hp) 



(38) 



(39) 



where kp (hp) is the Fermi wavenumbcr for the heavy 
(light) holes. This intrinsic spin Hall effect is driven 
by the band structure and its Berry phase connection, 
in sharp contrast to the extrinsic mechanisms previously 
proposed by Dyakonov and Perel in 1971 [7l| and oth- 
ers [72|, [73[. Sinova et al. independently proposed the 
intrinsi c sp in Hall effect for n-GaAs using the Rashba 
model 74j. Note that eq. (|3"8")) looks natural when one 



remember that the SU(2) gauge potential Af — Ei a iEt is 
coupled to the spin current as discussed in section I. How- 
ever, the topological structure of the Bloch wavefunction 
is needed to realize the intrinsic spin Hall effect. 

There are several experimental reports on the spin Hall 
effect. In semiconductors, the spin accumulation near the 
edge of the sample is measured to confirm the spin Hall 
effect. Kato et al. [75| observed the spin accumulation 
at the edge of the n-type GaAs in terms of the Kerr 
rotation. Wunderlich et al. (ZfiJ observed the circular po- 
larized LED from the recombination process of the holes 
in the p-type GaAs. Both experiments seems to be rele- 
vant to the spin Hall effect, but its origin, i.e., intrinsic 
or extrinsic, is still controversial ffl\ . The experiment 
on metals usually measures the inverse spin Hall effect, 
i.e., the voltage induced by the injection of spin current 
from the ferromagnetic metals [78[ . For some of the met- 
als, the first-principles band structure calculations are 
applied to predict the intrinsic spin Hall effect in e.g. Pt, 



which agrees reasonably well with the experimental ob- 
servations [7^]. The readers are referred to a review [77j 
for more details. 



III. OTHER INTERESTING SYSTEMS AND 
CONCLUSIONS 

There are several important and intriguing systems 
from the viewpoint of EEMF which arc not covered in 
section II. One natural question to ask is "What is the 
global aspect of the sub-Hilbert space and gauge fields, 
and their implications on the physical properties ?" . Of 
particular importance from this viewpoint is the topo- 
logical insulators [8(J. One can imagine the 2 copies of 
quantum Hall systems with opposite chiralities for up 
and down spins. This systems will show the quantized 
spin Hall conductance due to the nonzero Chern num- 
bers Cf = —C^. The spin Chern number is defined as 
C-\ — Cj,, which can be finite even for the T-symmetric 
systems. However, in the presence of the SOI, usually all 
the components of the spin are not conserved, and hence 
the spin current and Chern number are not well defined. 
Actually, the spin Hall conductance in spin Hall insula- 
tor, where the chemical potential is inside the band gap 
but still the spin Hall conductance is finite due to the 
band inversion, is not quantized (82T ? ]. Kane-Mele dis- 
covered a T-symmetric model with SOI, which shows a 
stable nontrivial topological phase with the helical edge 
channels even with the mixing of the spins [80j . This 
phase is characterized by the Z2 topological number in- 
stead of the Chern number. This idea is now generalized 
to 3-dimensions and also to arbitrary dimensions to clas- 
sify all the possible topologically nontrivial states as far 
as the band structure can be defined [83| . The basic con- 
cept is that the Bloch wavefucntions as the fiber bundle 
can have the globally nontrivial structure characterized 
by the index defined as the integral of the Berry curva- 
ture or EEMF over the 1st BZ. Combining the mapping 
between the different dimensions and symmetry class en- 
ables the classification of the topological classes including 
both the insulators and superconductors with gap. Re- 
cently the generalization to include the k- and r-spaces 
has been achieved also [lij . The readers are referred to 
recent reviews for more details [80j | . 

Noncentrosymmetric systems with SOI are also inter- 
esting laboratory to study the topological nature and 
EEMF. Rashba system is the representative one de- 
scribed by the Hamiltonian 



Hr = h Xe z ■ (s x p), 



(40) 



where e z is the direction of the potential gradient which 
breaks /-symmetry, s and p are the spin and momen- 
tum operators, respectively. This interaction has been 
experimentally demonstrated for the 2-dimensional elec- 
trons at the interface of GaAs system I8al or at the in- 
terface electrons in oxides system [86|, |83|, and also the 
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bulk 3-dimensional material BiTel [8S( | . The spin split- 
ting at each k-point means that the spin and the veloc- 
ity is tightly coupled, and hence the transport and mag- 
netic properties. Therefore, the various cross correlations 
between the magnetic and transport degrees of freedom 
are expected and actually observed experimentally in the 
Rashba systems, which constitute the important part of 
the spintronics. The spin Hall effect, spin Galvanic effect, 
and dynamical magneto-electric effect are the examples 
of this cross-correlation driven by the Rashba interac- 
tion [H, [9(|. Combined with the superconducting or- 
der parameter, the Rashba type SOI leads to even more 
interesting phenomena [9lJ. Because of the absence of 
the inversion symmetry, the classification of even and 
odd parity does not work there, and hence the singlet 
and triplet pairings are mixed. This fact leads to sev- 
eral unusual features such as the H c i beyond the Pauli 
limit [9l|. Furthermore, it is shown that the topologi- 
cal superconductor with helical Majorana edge channels 
as the Andreev bound states can be realized when the 
p-wave pair ing is dominant over the s-wave pairing com- 
ponent [92] ]. When the magnetic field or exchange field 
is applied to open the gap at the band crossing and the 
Fermi energy is within that gap, the spinless pairing su- 
perconductor is realized, and the Majorana fermions are 
expected to appear there [93j. These are consistent with 
the more general classification scheme of the topological 
superconductors [Hj]. 

In this article, we have discussed the physics of emer- 
gent electromagnetic field (EEMF) in condensed-matter 
systems. The gauge connections and fields are naturally 



introduced in various situations when the wavefunctions 
are constrained on some manifold in the Hilbert space. 
In the band structure and/or in the mean field theory, 
this gauge fields are not dynamical but frozen one or are 
controlled as parameters by the external electromagnetic 
field (EEM). However, the physical phenomena driven by 
EEMF are even richer than those of EMF because the 
lattice gauge theory, non-Abelian gauge fields, higher di- 
mensions, topological terms, are often relevant to EEMF. 
In this sense, the gauge fields can be a guiding and uni- 
fying principle in condensed-matter physics. Needless to 
say, the dynamical nature of EEMF in correlated sys- 
tems and spin liquids, which was not covered in this 
article, also remains to be an important subject, and 
continue to be a central issue in the condensed-matter 
theory 0, H, [T3l - 
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Fig. 1 Parallel transport of a vector on the sphere. Going around a loop, 
the direction of the vector has been changed corresponding to the solid angle 
enclosed by the loop. 
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Fig. 2 Scalar spin chirality as the gauge field for the conduction. When the spin 
directions S's form the solid angle Q, the conduction electrons whose spins 
(red arrows) are coupled to S's feel Q/2 as the gauge flux. 
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Fig. 3 Cluster model for the spin current mechanism of the electric polarization. 
The red arrow represent the direction of the flow and the orange arrow that of the 
Spin polarization of the spin current, which produces the electric polarization P 
through the spin-orbit interaction on the d-orbitals. 
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Fig.4 Experimental observation of the 
multiferroic behavior of DyMn03. (a) and 

(b) show the temperature dependence of 
the spontaneous electric polarization along 

c and a axes, respectively, for several values 
of magnetic field. There occurs the flop of 
The polarization under the magnetic field. 

(c) and (d) show the corresponding 
temperature dependence of the dielectric 
constant along a-axis. 
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Fig. 5 Anomalous Hall effect in SrRu03. 
(a)magnetization, (b) resistivity, 
(c) transverse resistivity as functions of 
temperature, (d) shows the transverse 
conductivity as a function of the spontaneous 
magnetization compared with the first 
-principles band structure calculations. 
The non-monotonous dependence is due 
to the magnetic monopoles in momentum 
space. 
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Fig. 6 Distribution of the Berry curvature in momentum space as a function 
of (kx,ky) at kz =0 calculated for SrRu03. The sharp peak at T-point is due to 
the band (anti)crossing acting as the magnetic monopole. 
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Fig. 7 Finite frequency anomalous 
Hall conductivity of SrRu03 in the 
Telahertz region. The observed 
behavior is well fitted by the simple 
model of the anti-crossing near the 
Fermi energy. 
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Fig. 8 Anomalous Hall conductivity as a function of the longitudinal conductivity 
calculated for a two-dimensional model with the anti-crossing near the Fermi 
energy ( pink curve). Various experimental data are also plotted. 
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Fig. 9 Anomalous Hall effect in NCI2M02O7. (a) The crystal structure contains the 
two inter-penetrating networks of tetrahedrons of Nd and Mo atoms, which are 
coupled antiferromagnetically. Due to the strong spin anisotropy Nd spins forms 
the non-coplanar configuration with scalar spin chirality at low temperature, which is 
transferred to conduction electron spins on Mo atoms as shown in (b). (c) The 
magnetic field dependence of the Hall resistivity for several temperature. Due to the 
scalar spin chirality, it shows the rapid increase as the temperature is lowered, (inset). 
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Fig. 10 Skyrmions and Skyrmion crystal in (Fe.Co)Si. Left panel: experimental 
observation by Lorentz microscope ((a),(b),(c)) at each point shown by the star 
In the phase diagram (g). Right panel: the corresponding figures obtained by the 
Monte Carlo simulation for DM magnet in 2D. 




Fig . 1 1 A single Skyrmion viewed as the mapping from the 2D real space to the 
unit sphere. This mapping is characterized by a topological number called Skyrmion 
number, i.e., how many times the mapping wrap the unit sphere. 




Fig. 12 Thermal Hall effect in Lu2V207. (a) the distribution of the DM vectors for a 
tetrahedron, (b) the magnetic field dependence of the thermal Hall conductivity at 
T=20K compared with the calculation (eq.(36) in the text). 



